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Abstract. In this paper a modified version of the Dynamic Renormalization Group 
(DRG) method is suggested in order to cope with inconsistent results obtained when 
applying it to a continuous family of one-dimensional models. The key observation is 
that the correct fixed-point dynamical system has to be identified during the analysis 
in order to account for all the relevant terms that are generated under renormalization. 
An application of this approach to the nonlocal Kardar-Parisi-Zhang equation resolves 
the known problems in one-dimension. Namely, obviously problematic predictions are 
eliminated and the existing exact analytic results are recovered. 
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Fluctuating surfaces appear in a wide variety of physical situations and have been of 
great interest in the last two decades [IllSlE]. These and other systems far from thermal 
equilibrium pose a major challenge in contemporary statistical physics. Behavior out- 
of-equilibrium is far richer than at equilibrium, and many intriguing scaling phenomena, 
such as self-organized criticality [4], or phase transitions between non-equilibrium 
stationary states [1], have been observed for long. However, despite the considerable 
achievements, the theoretical comprehension of non-equilibrium phenomena remains 
much poorer than our understanding of equilibrium phenomena. 

The renormalization group (RG), proven useful to explain universality in 
equilibrium continuous phase transitions, has also allowed some progress in 
understanding systems out-of-equilibrium. Nevertheless, in many cases the information 
RG analysis offers in not complete and limited to a certain range of dimensions. A 
classical example is the Kardar-Parisi-Zhang (KPZ) equation [3] where the Dynamic 
Renormalization Group (DRG) approach agrees with the analytic exact result in one 
dimension [1] but unable to provide results for the strong coupling phase in higher 
dimension. This clearly indicates that internal problems exist in the DRG calculation 
for d > 1. Actually, a remarkable result of Wiese [5] shows that the shortcoming of DRG 
in the KPZ system is not an artifact of a low order calculation (so called "one loop" 
calculation), but rather intrinsic to the method and extends to all orders. This situation 
motivated the development of other methods to deal with the KPZ system such as a 
scaling approach Self- Consistent Expansion (SCE) [7], Mode-Coupling [8] and others 
that were able to provide predictions for the exponents in more than one-dimension. 

A decade ago, a family of nonlocal growth models have been introduced in [9], 
known as the Nonlocal KPZ (NKPZ) equation, to account for nonlocal interactions in 
a system of deposited colloids, giving rise to roughness larger than the one predicted 
by the classical KPZ case. The authors studied the white noise case that was later 
generalized to spatially correlated noise in pilj. To be more specific, the equation they 
studied was 

dh(f,t) N f ,d ,^h(f,t)-Vh(f',t) , 

= uW^h (f, t) + ^ j d'r' ^^_^^d-p + V (r, t) , (1) 

where rj (r, t) is a noise-term modeling the fluctuation of the rate of deposition, which 
has a zero mean and is characterized by its second moment 

(r/(f,t)77(r^t)) = 2Do|r-r-^|"5(t-t') , (2) 

where d is the substrate dimension and Dq specifies the noise amplitude. Both papers 
[9l [To] investigated this problem using Dynamic Renormalization Group (DRG), and 
derived a complex phase diagram. Focusing on the strong coupling solution (in the 
KPZ sense [Ill3j) both papers found 

_ {d-2-2p) (d-2-3p) 
^-^^ (3 + 2-P)rf-6-9p ' 
where z is the dynamic exponent. The roughness exponent, a, characterizing the 
long distance spatial behavior, is obtained using the modified Galilean scaling relation 
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a + z = 2 — p. 

Unfortunately, the DRG results for the exponents summarized in Eq. ([3]) above, 
were found to be inconsistent with an exact result available in ID [12J predicting 
z = {3 — 3p)/2 when p = 2cr. A more systematic study using the Self- Consistent 
Expansion [llj led to the hypothesis that DRG fails in recovering the exact one- 
dimensional result since it does not account for new modes of relaxation generated 
by the special nonlinearity. This suggests going back the old Renormalization idea of 
identifying the right fixed point dynamical system around which the expansion should 
be is performed. The fixed point dynamical system is not necessarily of the same form 
as the original system, as is implicitly assumed by the standard DRG procedure. 

In this paper a modification of the standard DRG procedure that goes along those 
lines is suggested. This approach makes DRG more flexible, and succeeds in recovering 
the exact result for the case of NKPZ. Not less important, this approach could be useful 
in implementing DRG in other situations where long-range interactions are present, such 
as those appearing in the context of wetting of an amorphous solid by a liquid [131 [H] 
and in in-plane tensile crack propagation in a disordered medium flbi [16] . The main 
motivation here is to make the first step towards extending the range of applicability 
of DRG in a field that suffers anyway from a lack of analytical tools, in order to allow 
further progress in systems out-of-equilibrium. 

To understand the origin of the difficulty, consider the one loop DRG. The 
renormalization procedure is most succinctly described through the Fourier momentum 
q and frequency u modes, in terms of which Eq. ([T]) becomes 



h [k,uj\ = Go [k,LJ ] r] [k,uj 



h i k,uj 



(4) 



where Gq {k,ujj is the bare propagator given by Gq {^k,ij] = l/{uQk'^ — ito), and 
is a nonlinear functional of the height given by 



N 
N 



h yk, u 
h (k,uj 



1 



Gn { k,uj] k 




d'^qdVL 



jq- ik — q] h {q,il) h ik — q,iu — n) , (5) 



2 " V W J J (2vr) 
The one loop expression for the dressed propagator defined hj G [k,uj] = 
h ik, Lu) /t] ik,uj] is given [U [17] by 



G{k,uj) =Go{k,co)+A[-^) Gl{k,u 




d'^qdQ 



Y2Doq 



-2a 



k — q 



-p 



q- [k — q 



k- 



-q) ■ k 



xGo{k- q, uj-il^Go {q, n) Go (-g, -n) , (6) 
which, after some algebra (see appendix B in Ref. [Ij for example) gives 



G{k,0) =Goik,o) + ^^Gl[k,0) k^-PR, 



d + p-2-2a 
' Ad 



dqq'^ 



-3-20--P 



(7) 
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where = S^/ {27tY, and Sa is the surface area of a (i-dimensional unit sphere. 

In local KPZ the last equation is used to calculate the renormalization of the surface 
tension z/q. However, a look at equation highlights the problem. While the first term 
on the RHS scales as the second scales as k~'^~''. Here a distinction between three 
cases should be made: (a) When p < the correction term (proportional to A^) is 
irrelevant compared to the first term in the limit of small momentum (i.e. in the limit 
of large scales), (b) When p = both terms have the same scaling dimension. This 
situation is actually the case in the classical KPZ equation (with correlated noise), which 
is well studied for example in Refs. [171 [13 [IS]- And (c) when p > the correction 
is dominant over the first term. This means that in this situation the perturbative 
expansion produces more relevant terms than those originally present in the equation. 
More specifically a fractional Laplacian is produced under the renormalization. This 
implies that the fixed-point system in the space of dynamical systems is a-priori not the 
original model, and one needs to consider a more general form which contains such a 
term in the equation in the first place. Adding a vik^'^ term by hand and going through 
the same process, a new (partially) dressed propagator Gi {k,uj^ = 1/ {uik'^^P + iu) is 

obtained. Repeating the steps described above gives a 2"'^ order expansion for the full 
propagator, similar to Eq. ([7]) 



G {k,0) = Gi{k,0 



^Gl (K^^e-^K.'-^' 



2 + 2p 



Ad 



dqq' 



.d-3-2cr+p 



(8) 



This time, all the terms have the same scaling dimension so that the perturbative 
expansion is meaningful in the sense that higher order corrections are not more relevant 
than lower order ones. This allows to calculate the renormalization of the effective 
surface tension i>i when p > 0. 



Ui = Ui 



X'^Do d-2a-2 + 2p 



d-3-2a+p 



(9) 



^1 Ad 

Next, the renormalization of the noise term is calculated. The effective noise D is 
defined as the contraction of two terms according to 



'h l^k, ujj h l^k', u'j j = 2DG i^k, uj G [k', u'j 6'^ [k + k' ) 6 {u + u') .(10) 
The one loop expansion yields now 

d'^qdn 



2Dk 



-2a 



2Dok~^'' + 2 {2Do 

2 




(2rc) 



d+1 



X 



q- [k — q 



-2u 



-2a 



k-q Gx[k-q,uj-9.] \Gx{q,9L)\. 



(11) 



Notice that G\ \ k^ j was used in the last expression. In case p < this should be 
G^Vk^uj] as in standard DRG. Evaluating the integral in Eq. fill I) one obtains 



Dk 



-2a 



dqq' 



.d-3-4a+3p 



(12) 
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(or dqq'^^^^'^" when p < ). As before, the behavior of this equation is comphcated 
by the fc-dependence of the correction term, and there are three options: (a) When 
p < a the correction is irrelevant and the noise amphtude does not renormahze. (b) For 
p = a the second term is of the same order as the first term, and therefore renormahzes 
the noise amphtude (the KPZ equation is an example of this case since there p = a = 0). 
And (c) for p > a the correction term is more relevant than the first term. This means 
that in this situation the perturbative expansion produces an additional correlated noise 
which is more relevant than that originally present in the equation. This implies, just 
like for the propagator above, that the fixed-point system in the space of dynamical 
systems does not have the form of the original model, and a more general form with 
a new noise term Dik~^^ is considered. Doing the RG calculation from the beginning 
gives 



1 + ^ ■ ^ / dqq 



(13) 



uf 4 

Last, the one-loop contribution to the vertex Ap is calculated. Without getting 
into all the details, the final results is that the vertex does not renormalize to one-loop 
order Ap = Ap, since the structure of the perturbation theory is analytical in nature 
and cannot generate singular terms that renormahzes Ap. This is simpler than the 
non renormalization of the vertex in the classical KPZ case (with p = 0), where the 
correction is identically zero because of some exact cancelation of terms (see Fig. B.3(c) 

in m 

Following the standard rescaling procedure the following flow equations are 
obtained, 

dun 



dui ( ^\'^Dod-2a-2 + 2p\ , 



dn 







di 



Do{z -2a- d + 2a) p < a 



and 



-^ = D,[z-2a-d^2p^^-^] p>a, (15) 
dX 



Xp{a + z-2 + p). (16) 

The last step is a discussion of the complete RG flow for the NKPZ equation. Four 
sectors in the p, cx-plane, in which solutions can be looked for, are identifled. In the 
following, a detailed analysis in one of the sectors is presented, and results for the other 
sectors are provided. Sector I is deflned by p > and p < a. In this sector the flow 
equations are (fT^ b. (fT5!) a and (fT6|) . As traditionally done (in Refs. [HIH] for example), 
it is simpler to combine the flow equation into one equation for the coupling constant 
deflned here a,s g = K^X'^pDo/ufd. The RG flow of g becomes 

^^ = {2-d-p + 2a)g + 3{d-2a-2 + 2p)g\ (17) 
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and so the Fixed Points (FP) for g are 

^° = ° ^ = 3(2-^-2p + 2a) - ^'^^ 

A special dimension comes out of the last expression, the so-called critical dimension 
which is dc = 2 — p + 2a. In Fig. [1] the RG flow of the coupling constant g for various 
dimensions is presented. 




Figure 1. Coupling constant flow for the NKPZ equation in Sector I (p > and 
p < a), where dc = 2 — p + 2a. The three cases (a)-(c) cover all possible dimensions. 



(a) The case d < dc — p: As can be seen in Fig. ^a), the nontrivial FP is the only 
stable FP in this region, and by plugging it into the flow equations gives the following 
scaling exponents a = (2 — (i— p+2cr)/3 and z = ((i+4 — 2p — 2(t)/3. These exponents are 
the generalizations of the exponents of the classical KPZ system with spatially correlated 
noise [HI [IS [19] . 

(b) For dc — p < d < dc (Fig. [11(b)), the trivial FP is the only possible FP in 
the physical range as the nontrivial FP is negative. However, g^ is unstable, and so the 
system flows towards g = oo. However, just like the strong coupling regime in KPZ, it 
is inaccessible to a perturbative consideration, and one can just indicate its existence 
without having a quantitative prediction for its scaling exponents. 

(c) Last, as can be seen in Fig. [11(c), for d > dc the trivial fixed point (7q = is 
stable, and the exponents are given bya = {2 — d — p + 2a) /2 and z = 2 — p. These 
exponents correspond to the exponents of the Fractal Edwards- Wilkinson equation (i.e. 
a linear equation with a fractional Laplacian) with correlated noise. In addition, for a 
higher bare value of the coupling constant g > g* , the system flows to (7 = cxd, which 
signals the appearance of the strong-coupling regime, again inaccessible to a perturbative 
approach. Thus, in this range of dimensions there is a possible phase transition between 
a weak-coupling to a strong-coupling regime. 

The analysis of the possible phases in the other three sectors follows the same lines 
and the results are summarized in Table [H Note that for the strong coupling phases, 
no quantitative result for the corresponding exponents is possible, apart from pointing 
out their existence. 

As can be appreciated, the full description of the results given in Table [H is quite 
rich. In order to gain more insight into this a special attention is given to the interesting 
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Table 1. A complete description of all the possible phases of the NKPZ problem using 
the modified DRG scheme for any value of d, p and a. The first two columns give the 
scaling exponents a and z for a particular phase, and the third column states each 
phase's validity condition. The values of the scaling exponents in the strong coupling 
phases are not accessible, and one can just indicate their existence. It is not even 
known whether they correspond to the same phase, and thus described by the same 
exponents or not. 

a z validity 

p>0, p<a,d>2-p + 2<7 
p>0,p<a,d<2-2p + 2(T 

p>0, p<a,d>2-2p + 2a 

P > 0, p > (7, d < 3/2 

P > 0, p > (7, d > 3/2 

p>0, p>(T,d>2 + p 

p < 0, p < a 

p<0, p>a,d>2 

p < 0, p > a 



{2-d- p + 2(t)/2 
{2-d- p + 2(j)/i 
Strong Coupling - sector 1 

(2-d)(2-d+p) 
2(3-2d) 

Strong Coupling - sector 11 
{2-d- p + 2(7)/2 
{2-d + 2(j)/2 
{2-d)/ 2 

Strong Coupling - sector IV 



2- p (FCEW) 

(d + 4-2p-2CT)/3 



9 _ „ _ (2-d)(2-d+p) 
^ 2(3-2d) 



2- p (FCEW) 
2 (CEW) 
2 (EW) 



one dimensional case with a > 0, where the following dynamic exponent is obtained: 

? w '"^ ■ (19) 

(3 - 3p)/2 p > ^ ' 

Note that for p > the exact one-dimensional result is recovered [12] , as was suggested 
in [llj, and thus a major problem with DRG described above is solved. 

To summarize, in this paper a modification of the classical DRG approach to 
systems out of equilibrium is presented in order to resolve problems with the results 
derived using traditional DRG for the Nonlocal KPZ equation. This approach extends 
beyond the NKPZ system, to any system out-of-equilibrium that produces under 
renormalization relevant terms which are not present in the original model. For the 
NKPZ system (as well as for the Nonlocal Molecular Beam equation [20] and the Fractal 
KPZ equation [21j) it is found that for certain values of the parameters a fractional 
Laplacian is generated under renormalization (p > 0), or a correlated noise term {p > a). 
Thus, an inclusion of these terms in the original model (or put differently, by considering 
the right fixed-point dynamical system) leads to a correct description, and resolves the 
above-mentioned inconsistency with the exact ID result when p = 2a. 

The important lesson from this discussion is that in order to obtain a correct 
description of a given problem using DRG, one needs to verify that the correct fixed 
point dynamical system has been identified, and perform the perturbative expansion 
using this fixed-point system instead of the original model. Following this lesson can 
help to extend the range of applicability of DRG, since the direct contradiction with the 
exact result is settled. 

An interesting application of these ideas could be the case of a driven wetting 
line of a fluid on a rough surface [131 E] o'^ the mathematically similar problem of an 
in-plane tensile crack propagating in a disordered material [IHl HE] (a moving rather 
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than a pinned interface). In these problems, a long-range interaction term exists at the 
nonlinear order, and it is therefore vulnerable to similar difficulties. As suggested in 
Refs. [m [16], it could be that that the Edwards- Wilkinson system [2] is the relevant 
fixed-point system for the rough phase of these physical models. More work in that 
direction is needed to clarify this issue. 
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